Dewetting of solid films with substrate mediated evaporation 
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The dewetting dynamics of an ultrathin film is studied in tiie presence of evaporation -or reaction- 
of adatoms on the substrate. KMC simulations are in good agreement with an analytical model with 
diffusion, rim facetting, and substrate sublimation. As sublimation is increased, we find a transition 
from the usual dewetting regime where the front slows down with time, to a sublimation-controlled 
regime where the front velocity is approximately constant. The rim width exhibits an unexpected 
non-monotonous behavior, with a maximum in time. 



I. INTRODUCTION 

Dewetting, the process by which a thin film breaks 
up into droplets, is a non-equilibrium process driven by 
the reduction of the total energy. A large number of 
recent experimental studies have been devoted to the 
analysis of the dewetting of ultra-thin solid films [iHl, 0- 
[^. The main mass transport process during dewetting 
is surface diffusion. However, early work from Srolovitz 
and SafranQ on solid dewetting have already pointed out 
the role of sublimation, which cannot be avoided at high 
temperatures. These authors have discussed the case of 
sublimation directly from the film to the vapor. Here, we 
discuss the regime where atoms from the film diffuse on 
the substrate, and may evaporate when they arc on the 
substrate. 

Further motivation of the present study comes from 
recent work on the dewetting of SOI systems (Si/Si02), 
which have revealed the presence of reaction of Si on the 
substrate (with subsequent sublimation of the products 
formed in the reaction) ■ Such a phenomenon could 

be modeled by a simple sublimation rate to a first approx- 
imation, neglecting the substrate shape changes induced 
by the etching related to the reaction. In addition, subli- 
mation could also be relevant for the dewetting of metal- 
lic films where the temperature can be raised close to the 
melting temperaturefJQ- 

In this article, we analyze the detailed dynamics of 
dewetting with the possibility for atoms to leave the film, 
diffuse on the substrate, and evaporate. We perform Ki- 
netic Monte Carlo (KMC) simulations, and provide a 
model inspired from that of Ref.[ll|, which is in good 
agreement with KMC simulation results. We show that 
substrate evaporation leads to a novel dynamical regime 
where the film edge velocity is approximately constant. 
In addition, the rim width exhibits an unexpected non- 
monotonous behavior in time, reaching a maximum value 
at some intermediate stage. 



II. KMC SIMULATIONS 

We model the dewetting of a crystalline film using a 
solid-on-solid model on a 2D square lattice. The sub- 



strate is fiat and frozen, and is associated to height z = 0, 
whereas the film has local height z > 1. The local height 
will change as atoms diffuse or evaporate. There is no 
incoming ffux of atoms in this model. 

An atom at the surface can hop to nearest neighbor 
sites with rates z^„, when it is not in direct contact with 
the substrate, and r„ when it is in contact with the sub- 
strate (z = 1). In our model an atom has to break all its 
bonds to hop. The hopping barrier is therefore given by 
the binding energy of the atom. An isolated atom over 
the substrate can also evaporate with rate r^.. Then: 
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where v is an attempt frequency, T is the temperature 
(in units with = 1), n'm the number of in-plane near- 
est neighbors of the atom before the hop, J is the lateral 
bond energy between two atoms in the film, Jq is the 
vertical bond energy between two atoms in the film, E^g 
is the energy barrier for desorption of an atom in contact 
with the substrate and Eg is the adsorbate-substrate ex- 
cess energy. The model is presented in Figl. We choose J 
as the energy unit, so that J = 1 in the following. Since 
Jo appears in all rates in Eas. (|ll2l3p . we choose v^^ as 
the time unit, where 
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Hence, in the following, we set i^q ~ 1. 

The parameter £^5 controls the wetting properties of 
the film on the substrate. For example, when Es < 0, 
the film completely wets the substrate. In this work, we 
consider Es > (partial wetting). It can be shown [l2| 
that the surface excess energy Es can be written as 



Es — Eav + Eas — E, 



sv, 



(5) 



where Eav i Eas a-nd Esv a-re the energy per site of the 
adsorbate-vacuum surface, adsorbate-substrate interface 
and substrate-vacuum surface, respectively. This param- 
eter is then identical (with opposite sign) to the spread- 
ing coefficient S defined in references Ref. [ij, [l^ for the 
study of liquid dewetting. 
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FIG. 1: (Color online) Schematics of the Kinetic Monte Carlo 
(KMC) model, with hopping rates Vn and r^, and sublimation 
rate rg. 



The algorithm used in the simulations is the following. 
Wc list all atoms into classes. Each class is characterized 
by the number of in-plane neighbors n of the atom and 
by the existence or not of a nearest neighbor belonging 
to the substrate. At a given time t, we calculate the 
probabilities per unit time Wi of all possible events (an 
event is the motion of an atom originally at position i) 
given either by Eq.©, Eq.© or Eq.©, and the sum W 
of all those rates (for all mobile atoms). We increment 
the time by a 6t, which is equal to the inverse of the sum 
of the rates of all possible events [l^, \/W. We choose 
the event with probability Wi/W . In the case of diffusion, 
the atom moves with equal probability in any of the four 
possible directions. 

The initial configuration is a film of constant height 
covering the substrate, leaving only a small stripe uncov- 
ered along the (10) direction. The film edges are straight 
in the beginning. Atoms diffuse at the surface of the 
film and on the substrate. However, they can evaporate 
only when they are on the substrate. If we let the system 
evolve for a long enough time, all atoms should evaporate, 
and there will be no atoms left on the substrate in the 
final state. However, we are only interested in the early 
stages of the dynamics, where the dewetting dynamics 
occurs. As shown in Figl2l the dewetting dynamics pro- 
ceeds with the formation of receding dewetting rims at 
the edges of the film. 

We varied the sublimation energy barrier E^s to in- 
vestigate the influence of sublimation in the dewetting 
regimes. If E.^s/T is sufficiently high, we recover the 
case considered before [111, [l^l with no sublimation (the 
adsorbate dewets driven only by diffusion and attach- 
ment/detachment of adatoms). 



III. MODEL 

A. Rim width and position dynamics 

In order to analyze the dewetting dynamics observed 
in KMC simulations, we use a ID model in which the 
dewetting front is assumed to be straight, and invariant 
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FIG. 2: (Color online) Snapshots of the dewetting dynamics 
with substrate evaporation observed in KMC. (a) We start 
from an initial trench in the film. (b,c,d) Two dewetting rims 
form and move in opposite directions. Note that the rim 
width exhibits a maximum in the intermediate stages in (c). 
Simulation parameters: Evs = 5, Lx = 400, Ly = 800, T — 
0.4, Es = 0.5 and h — 3. The z height scale of the surface is 
indicated on the right. 



along the direction y. The position of the edge of the 
film is denoted as xi. At xi, the height exhibits a jump 
of height hi between the substrate and the rim top facet. 
On the other side of the dewetting rim, a bunch of atomic 
steps is at position X2 > xi. On the film for x > X2, 
the height stays at the value of the initial film height h. 
Defining /i2 = hi — h, we have: 

hidfXi = —^Ddxc\i + {x2 — xi)dthi — i}DsdxCs\i(^6) 
h2dfX2 = -rtDdxc\2, (7) 

where c and D are the adatom concentration and diffu- 
sion constant on the rim top facet, and Cs and Dg are 
the adatom concentration and diffusion constant on the 
substrate. In addition, fl is the atomic volume. 

On the substrate, the adatom concentration Cg obeys 
a quasistatic diffusion equation with evaporation 



Cs 
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where Dg and Tg are the diffusion constant and subli- 
mation time of adatoms on the substrate. The general 
solution of Eq.® reads: 



Cs = ^cosh(a;/a;s) + B s.mh{x / s s) , 



(9) 



where Xs = {DgTsY^'^ is the sublimation length. We 
consider a geometry similar to that of KMC simula- 
tions: Two dewetting fronts, separated by the distance 
is = 2x1, recede in opposite directions. The boundary 
condition at the fronts is assumed to be instantaneous 
attachment-detachment kinetics, leading to: 



Cs 1 



Using this condition in Eq.Q, we find 

J, cosh(a;/a;s) 
"^^«cosh(4/2x,)' 



(10) 
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and, as a consequence the contribution of evaporation to 
the dewetting front velocity in Eq.(l6l) is: 
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Vs = nDsd:cCs\i = n y—j <qtanh(a;i/a;,,). (12) 

On the top of the rim, the adatom concentration c 
obeys a quasistatic diffusion equation 



Ac = 0. 



(13) 



Assuming translational invariance in the other direc- 
tion y, and instantaneous attachment-detachment kinet- 
ics leading to fixed concentrations ci = Cgq exp{Es /hiT) 
and C2 = Ceq at xi and X2, we obtain for xi < x < X2- 



ci 



X — Xl 

X2 - a;i 



(c2 - ci), 



(14) 



so that dxc\i — dxc\2 ~ (c2 — ci)/{x2 — xi). The dynam- 
ical equations may then be re-written as: 

no C2 - ci dthi Vs 

Otxi = h(a;2-a;i)— ^ — , (15) 

III X2 — xi hi hi 



dtX2 = - 



Q.D C2 — ci 



/l2 X2 - Xi 



(16) 



Now hi and h2 = hi — h are time dependent. Moreover, 
C2 — ci also depends on time via hi. 



B. Rim height evolution 

Since we assume no evaporation on the film, we use the 
same theory as in Rcf. [11] for the rim height evolution. 
Following Ref. [l3| we assume that the rim height evolu- 
tion is controlled by 2D nucleation on the facet with a 
nucleation rate: 
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where 7 is the atomic step free energy on the film. The 
monolayer islands formed by the 2D nucleation process 
then spread on the rim top facet with the velocity [ll|: 



(18) 



where Czip ~ 0.25 is a number. We will also assume 
that the temperature is high enough so that steps are 
isotropic, and 7 w 7. 

Assuming a total rim length L, the rim height evolu- 
tion reads! 



dthi = J min 



Zip 
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(19) 



which accounts for both regimes of single and multiple 
nucleation. 

The model equations (|15I16I19P are similar to that of 
Rcf. [m . The only difference is the term Vs which appears 
in Eq.l^. 



IV. COMPARISON BETWEEN MODEL AND 
KMC 

A. KMC model parameters 

In this section, we shall evaluate the model parameters 
as a function of the KMC parameters. 

On the film, the adatom concentration is given by de- 
tailed balance |ll|: 



Ceq ^ C ^ , 



and the diffusion constant reads: 



(20) 



(21) 



The energy cost for an adatom to be in contact with the 
substrate is Es, and therefore: 



n 



Es/T 



(22) 

(23) 



Finally, the sublimation time on the substrate reads: 

= e(^--^^)/^, (24) 
where Et,s is the substrate dcsorption energy. 

B. Quantitative comparison 

Wc have measured Xi, hi, and ^ as a function of time 
for various parameters in KMC simulations. As opposed 
to the case without sublimation, we have not found any 




t 

FIG. 3: (Color online) Front position xi for = 3 (red 
circles), 4 (green squares), 5 (blue squares), and 6 (magenta 
triangles). Other parameters: = 0.5, — 400, Ly = 400, 
T = 0.4, and h = 3. Symbols are results from KMC simula- 
tions, and the solid lines represent the numerical solution of 
the model of section Hill 
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FIG. 4: (Color online) Rim width for Evs = 3 (red circles), 4 
(green squares), 4.5 (dark green crosses), 5 (blue diamonds), 
5.5 (purple x's), and 6 (magenta triangles). Other parame- 
ters: Es = 0.5, = 800, Ly = 400, T = 0.4, and h = 3. 
Symbols are results from KMC simulations, and the solid lines 

□ 



(even partial) analytical solution of the model. We there- 
fore resort to a full numerical solution of Eqs. (jl5ll6ll9p 
The agreement between KMC and the numerical solu- 
tion of the model was checked for xi , i, and hi . As shown 
in Fig |3l4l5l the agreement is good (note that there is no 
fitting parameter in the model). Interestingly, the evolu- 
tion of the rim width is non-monotonous, and exhibits a 
maximum in time both in the model and in KMC simu- 
lations. However, the early-time dynamics is difficult to 
reproduce since the first layer nucleation is not described 
correctly within our model, as seen in Fig|31 In addition, 
the measurement of the rim width in the early stages of 
the dewetting process in KMC simulations is delicate. In- 
deed, we simply measure the area of the film with height 
z > h + 1, and divide by the total rim length 2Ly. Since 
adatoms and small monolayer islands on the film have 
z = h + 1, they are also counted. Hence, this procedure 
is inaccurate in the initial stages when the total rim area 
is small. 



V. DISCUSSION 
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FIG. 5: (Color online) Rim height for Evs ~ 3 (red circles), 4 
(green squares), 5 (blue diamonds), and 6 (magenta triangles). 
Other parameters: E^ = 0.5, = 400, Ly = 400, T = 0.4, 
and h = 3. Symbols are results from KMC simulations, and 
the solid lines represent the numerical solution of the model 
of section IIIII 



We have studied the dewetting dynamics of thin solid 
films in the presence of substrate sublimation. The over- 
all agreement between KMC and the analytical model 
indicates that we have caught the main physical ingredi- 
ents of the dynamics. 

A. Transition time 

In this section, we discuss the transition between the 
early time diffusion-limited regime, and the late time 
sublimation- limited regime. 



Since the nucleation rate J decreases exponentially 
with hi, the evolution of hi essentially occurs in a short 
transient regime at short times. We shall therefore base 
our analysis on the fact that hi evolves very slowly. 

In the competition between diffusion and sublimation, 
the process which controls the dynamics is the one which 
leads to the largest contribution to the front velocity. For 
fixed hi, diffusion leads to xi ~ t^^^, and sublimation 
leads to xi ^ t. Hence, we expect that diffusion always 
wins at short times, while sublimation controls the late 
time regime. 

Let us first consider the diffusion limited regime (where 
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FIG. 6: (Color online) Rim width given by the model (solid 
line) for E^s = 5, Es = 0.5, L^^ = Ly ^ 400, T = 0.4, and 
ft = 3. The dashed line corresponds to loiff Ea. (|25|l . and the 
dotted line to IsuU Eq.JgeJ. 



Vs is negligible), with fixed hi, leading to 



[2Dn{ci - C2) {h-^ - h^^) t] ^''^ . (25) 



hi ^ h 



In the late-time dynamics, which is dominated by sub- 
limation, the front moves approximately at constant ve- 
locity. Hence, we expect the dewetting rim to be roughly 
stationary, with a slow variation of £. An evolution equa- 
tion for i is obtained by taking the difference between 
Eq. ([T5l) and Eg. ([T6l) . and assuming that dthi <^ Vg/i- 
Plugging dt£ — into this equation leads to 



iSub 



on 



Cl - C2 



hi — h 



(26) 



B. Sublimation from the film 

The situation discussed in the present work can be 
compared with the results of Srolovitz and SafranQ, 
which also find a front evolution with constant veloc- 
ity. In their model, there is simultaneously sublimation 
and adsorption between the film and the vapor, so that 
the net growth of a flat film vanishes. The motion of 
the film height h is then driven by local chemical po- 
tential variations. Since the local chemical potential is 
proportional to the surface curvature, we obtain in the 
small slope approximation: dth = Adxxh. Typically, we 
expect A = kpeqfl'js/kBT, where k is the evaporation- 
condensation kinetic coefficient, p^q is the equilibrium 
vapor density, and 75 is the surface stiffness (which is 
assumed to be isotropic). 

The edge of the film, which exhibits a larger curvature, 
has a larger chemical potential, so that the sublimation 
at the edge is reinforced and is not compensated by ad- 
sorption from the vapor. Therefore, the front recedes by 
sublimation at the film edge. Srolovitz and Safran find 
a front velocity V = Atand/h, where h is initial film 
height, and 9 is the contact angle, which is related to Es 
viafl^: 



Es/J^ ^{l^cos{e)). 



(29) 



As a summary, the modeling of Srolovitz and Safran cor- 
respond to the case of a film in equilibrium with its vapor, 
a situation which is very different from the one consid- 
ered in the present study. However, the dewetting pro- 
cess qualitatively obeys the same behavior, with fronts 
moving at constant velocitv(20|. 



Recalling that hi is not fixed, but in fact increases slowly, 
the expression of isuU should be interpreted as a slow 
decrease of i caused by the increase of hi. As shown on 
Fig|4l ^Diff and tsuU provide reasonably good approxi- 
mations of £ at short and long times. 

Combining the two expressions, we obtain the 
crossover time from the relation ^oiff = ^Subi as: 



Dnc,„ Es 



2vi 



T hi~h' 



(27) 



where we have assumed that Es ^ Thi. In addition, 
since the growth of hi slows down exponentially with 
time, we assume that h/{hi — h) is of order of 1. Finally, 
we obtain 



DD.Ceq E^ 



2vl 



s 
T ' 



(28) 



As a consequence, we see that the crossover time will 
decrease as sublimation (and hence Vs) increases. This 
formula provides the correct order of magnitude as com- 
pared to the maximum observed in Fig]?] 



C. Discussion of experiments 

Dewetting of SOI systems, consisting of Si crys- 
talline films on amorphous Si02 substrates, has been 
investigated recently The presence of reaction 

Si-t-Si02 — i-2SiO has been suggested in several experi- 
ments, and there is still an issue about the determination 
of the place where the reaction occurs. If Si atoms dif- 
fuse through the triple line on Si02, and then react with 
Si02, then the situation would be quite similar to our 
model system. (Note that the reaction should lead to an 
etching of the substrate, which is not taken into account 
in our model). This scenario is consistent with that pro- 
posed for the deoxydization of Si02 films on Si substrates 
jlTl [isj. However, the reaction mechanism is controver- 
sial, and a more complex scenario involving the diffusion 
of oxygen at the interface was also proposed[l3|. There- 
fore, we expect substrate reaction to lead to a change in 
the time evolution of the position of straight fronts, from 
a power law behavior with an effective exponent close to 
1/2, to an approximately linear behavior. This crossover 
could appear with increasing temperature if reaction is 
activated at high temperatures. 
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Substrate evaporation could also be important for the 
dewetting dynamics of metallic films where the tem- 
perature is often increased up to temperatures where sub- 
limation cannot be neglected. 

VI. CONCLUSION 

In conclusion, in the presence of adatom sublimation 
on the substrate, the motion of a dewetting front exhibits 
a crossover from a constant thickening and widening of 



the rim, accompanied with a slowing down (with the front 
position obeying an approximate power-law t^^^), to a 
regime where the velocity is approximately constant. In 
this latter regime, the rim height increases slowly, and 
the rim width decreases slowly. Surprisingly, the rim 
width exhibits a maximum at the crossover between the 
two regimes. We have obtained quantitative agreement 
between KMC simulations and a model with diffusion 
limited mass transport on the rim, nucleation limited in- 
crease of the rim height, and evaporation on the sub- 
strate. 
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